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The geometric aspects of quantum mechanics 
are underlined most prominently by the concept 
of geometric phases, which are acquired whenever 
a quantum system evolves along a closed path 
in Hilbert space. The geometric phase is deter- 
mined only by the shape of this path^ ^ and is 
— in its simplest form — a real number. How- 
ever, if the system contains degenerate energy 
levels, matrix-valued geometric phases, termed 
non-abelian holonomies, can emergd^. They play 
an important role for the creation of synthetic 
gauge fields in cold atomic gases^ and the descrip- 
tion of non-abelian anyon statistics^. Moreover, 
it has been proposed to exploit non-abelian holo- 
nomic gates for robust quantum computation^^^ 



OO In contrast to abelian geometric phases^^, non- 
^"^ abelian ones have been observed only in nuclear 
quadrupole resonance experiments with a large 
number of spins and without fully characterizing 
the g eome tric process and its non- commutative 
nature^^^. Here, we realize non-abelian holo- 
nomic quantum operations'^ '^ on a single super- 
conducting artificial three-level atom'^ by apply- 
ing a well controlled two-tone microwave drive. 
Using quantum process tomography, we deter- 
mine fidelities of the resulting non-commuting 
gates exceeding 95%. We show that a sequence of 
two paths in Hilbert space traversed in different 
order yields inequivalent transformations, which 
is an evidence for the non-abelian character of the 
implemented holonomic quantum gates. In com- 
bination with two-qubit operations, they form a 
universal set of gates for holonomic quantum com- 
putation. 

A cyclic evolution of a non-degenerate quantum system 
is in general accompanied by a phase change of its wave 
function. The acquired abelian phase can be divided into 
two parts: The dynamical phase which is proportional to 
the evolution time and the energy of the system, and the 
geometric phase which depends only on the path of the 
system in Hilbert space. This characteristic feature leads 
to a resilience of the geometric phase to certain fluctua- 
tions during the evolutioiP^Hi^, a property which has at- 
tracted particular attention in the field of quantum infor- 
mation processings^. However, universal quantum com- 
putation cannot be based on simple phase gates, which 
modify only the relative phase of a superposition state, 
unless they act on specific basis states^. Furthermore, 
geometric operations acting on degenerate subspaces have 



Oh 



> 

00 

in 

^' 

o 
m 



'i 



been proposed for holonomic quantum computation fully 
based on geometric concepts^. In this scheme, quantum 
bits are encoded in a doubly degenerate eigenspace of the 
system hamiltonian h{\). The parameters A are varied 
to induce a cyclic evolution of the system. When the sys- 
tem returns back to its initial state, it can acquire not 
only a simple geometric phase factor, but also undergoes 
a path-dependent unitary transformation, a non-abelian 
holonomy, which causes a transition between the eigen- 
states in the degenerate subspace. Because these trans- 
formations depend only on the geometric properties of 
the path, they share the noise resilience of the geometric 
phase. 

In the original proposal^, the parameters A are changed 
adiabatically in time to guarantee the persistence of the 
degeneracy. Adiabatic holonomic gates have been pro- 
posed for trapped ions^, superconducting qubits^^ ^^ and 
semiconductor quantum dots^^. However, they are dif- 
ficult to realize in experiment because of the long evo- 
lution time needed to fulfill the adiabatic condition. In- 
stead, Sjoqvist et al.l^ have proposed a scheme based on 
non-adiabatic non-abelian holonomies^^ which combines 
universality and speed and can thus be implemented in 
experiments. 

The main idea is to generate a non-adiabatic and cyclic 
state evolution in a three-level system which results in a 
purely geometric operation on the degenerate subspace 
spanned by the computational basis states |0) and |1). 
The third state |e) acts as an auxiliary state and remains 
unpopulated after the gate operation. This is achieved by 
driving the system with two resonant microwave pulses 
(Fig. Ilk) with identical time-dependent envelope ^^(t), 
but different amplitudes a and h satisfying \a\^ + \h\^ = 
1 (see Fig. Ilfc). The hamiltonian of the system in the 
interaction picture is 



h{t) 
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(a|e)(0|+6|e)(l|+/i.c.), 
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and it causes the initial basis vectors to evolve to the 
states \^IJi{t)) = Texp(-^Jh{t')dt'j\i) {ij = 0,1), 

where T denotes time ordering. In contrast to adiabatic 
schemes, the \ipi{t)) are not instantaneous eigenstates of 
h{t). By keeping the complex amplitude ratio a/b of the 
pulses constant, no transitions between states are induced 
and the evolution satisfies the parallel transport condi- 
tion, {ilji{t)\h{t)\iljj{t)) =0. As a consequence, the evolu- 
tion is purely geometric with no dynamic contributions. 

r 

If the pulse length r is chosen such that J Q{t)dt = 27r, 



the degenerate subspace undergoes a cyclic evolution, and 
the final operator which acts on the |0) and |1) basis states 
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FIG. 1. Geometric gate operation on a three- level sys- 
tem, a, Two drive tones with amphtudes a and b and pulse 
envelope Q{t) couple the |0) and |1) to the |e) state. To low- 
est order, the direct transition between |0) and |1) states is 
forbidden, b, Pulse sequence for process tomography: The 
input state is prepared by sequentially applying pulses on the 
|0) o |e) and |e) o |1) transitions (see Methods). The holo- 
nomic gate is formed by the simultaneous application of two 
pulses with envelope Q(t) and amplitudes a and b. A set of 
pulses on the |0) ^ |e) and |e) ^ |1) transitions followed by 
the measurement completes the process tomography, c, Holo- 
nomic gate represented on a fibre bundle. A path in the fiber 
bundle connects initial and final states. The projection of the 
path in the fiber bundle onto the base manifold spanned by 
the states {|'0i)('0i|}, forms a closed loop C along which the 
system is parallel transported. The difference between the 
initial and final points lying on a single fiber corresponds to 
the matrix valued holonomy U{C), an element of the unitary 
group U{N). 
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where we have parameterized the drive amplitudes by the 
relation &^ tan6>/2 = a/h. In a geometric picture, the dy- 
namics of the system can be visualized on a fibre bundle 
in which the closed loop C in the base manifold deter- 
mines the holonomic operation U{C) (Fig. lib). Different 
values of and (j) correspond to different paths C. 

In our experiments, we realize the holonomic gates us- 
ing a three-level superconducting artificial atom of the 
transmon type embedded in a 3D cavit}'^^ (Fig.[2p). 
The cavity is made of aluminium with inner dimensions 
32 X 15.5 X 5mm^. The frequency of the fundamental 
mode is uj^ol2i\ :^ 8.999 GHz as measured by transmis- 
sion spectroscopy using the circuit shown in Fig.[2]3. The 
transmon is made of two 500 x 250 \ivc? aluminium elec- 
trodes separated by 130 /im and connected via a Joseph- 
son junction (Fig.[2]3,d). It is oriented parallel to the 
electric field of the fundamental mode pointing along the 
smallest dimension of the cavity (Fig. [2k). The measured 




FIG. 2. Transmon qubit in a cavity resonator, a, Alu- 
minium cavity with an embedded sapphire chip containing 
two transmon qubit s. The left qubit is not used in the exper- 
iment. The electric field profile E of the fundamental mode is 
sketched in the upper part of the cavity, b, Microwave control 
pulses at the transition frequencies cjoe and cjei are created by 
modulating two continuous signals using in-phase/quadrature 
(I IQ) mixers. Modulation pulses in the / and Q channels are 
generated using arbitrary waveform generators. The control 
pulses and the readout pulse at frequency cjro are combined 
and fed into the cavity. The transmitted signal is amplified 
and detected in a heterodyne measurement at room tempera- 
ture and analyzed on a computer using a digitizer (ADC), c, 
Optical and d scanning electron micrograph of the employed 
transmon device. 



coupling strength between the transmon and the cavity 
field is ^/27r ~ 110 MHz. To read out the state of the 
transmon we measure the state-dependent transmission 
through the cavit}^. The transition frequencies mea- 
sured by Ramsey spectroscopy are a;oe/27r ^ 8.086 GHz 
and cjei/27r ::^ 7.776 GHz. The decay times of both ex- 
cited states are Ti = 7 ± 0.1 /is, and the dephasing times 
are T^"" = 8.0 ± 0.1 fis and Tf = 3.9 ± 0.1 /j.s. 

Different holonomic gates are realized by adjusting the 
ratio of the two drives a/b — e^^id.iiO/2 to values be- 
tween and 1/a/2 with (/) = 0. In geometric terms, this 
corresponds to a gradual deformation of the loop C in the 
projective Hilbert space resulting in a continuous change 
from the phase-flip gate a^ to the NOT gate cFx- The 
envelopes VL{t) are truncated Gaussian pulseP^ with a 
width of cr = 10 ns and a total pulse length of 40 ns. 
With pulses of this length, off-resonant driving of neigh- 
boring transitions is negligible. The performance of the 
holonomic gates is characterized by process tomography 
performed on all three levels (see Methods). The diago- 
nal elements of the reduced process matrix x ^i"^ shown 
in Fig. [3^ as a function of 6 with = tt. The experimen- 
tally obtained results are in good agreement with theory. 
For 6> = a single drive on the |e) ^ |1) transition causes 
a phase shift corresponding to xzz = 1 and the opera- 
tion U{0 = 0) = cr^. In the case = 7r/4, the Hadamard 




FIG. 3. Process tomography of holonomic gates, a, 

Diagonal elements xa of the process matrices for ideal (lines 
and experimental (dots) geometric gates as a function of 0. 
Black circles correspond to the trace of the reduced process 
matrix x- (t>) Bar chart of the real part of the reduced mea- 
sured process matrix Xexp of the geometric Hadamard gate H 
and (c) the NOT gate . The wire frames show the theoreti- 
cally expected values. The small (0.3%) imaginary parts of x 
are not shown. 



transformation H = (a^ — ax)/^/2 is generated (Fig.jSp). 
= 7r/2 corresponds to pulses with equal amplitude and 
generates a NOT gate ax (Fig. [3 3). Because of dephas- 
ing and relaxation of both excited states as well as finite 
fidelities of the microwave pulses, the |e) state is slightly 
populated after the gate operation. This leakage is quan- 
tified by computing the trace tr (x) ^ 0.96 of the reduced 
process matrix x (black dots). 

The fidelity of the Hadamard transformation is J^h ^ 
95.4% and the fidelity of the NOT operation is J- not ^ 
97.5%. The numerical solution of a master equation in- 
cluding dissipative processes results in fidelities of T = 
97.6% for both processes, in good agreement with the 
experimental values. 

To explicitly show that different loops of the state 
vector in the Hilbert space result in non-commuting 
quantum gates, we sequentially apply the geometric 
Hadamard and the NOT gate in alternating order. The 
non-abelian character of the operation yields either the 
operation NOT ■ H = -(icr^+ 1)/V2 (Fig.[4^) or H ■ 
NOT = {lay - 1)/V^ (Fig.Ho). We visualize the op- 
erations on the Bloch sphere (Fig. [IJ^-d): NOT • H cor- 
responds to a TT-rotation about the ilf-axis followed by 
a TT rotation about the X-axis. This is equivalent to a 
rotation about the F-axis by 7r/2. On the other hand, 
H ■ NOT corresponds to a rotation in the opposite direc- 
tion, which is in clear contrast to the former operation. 
By concatenation of operations other than the Hadamard 
and the NOT operations, rotations about arbitrary axes 
corresponding to a representation of the complete SU(2) 
group can be realized. 



FIG. 4. Non-commutativity of holonomic gates. Process 
matrices for NOTH (a) and HNOT (b) gates with fidelities 
of 95%. Because of the non-abelian character of the geometric 
operations the resulting processes are different. This can be 
visualized on the Bloch sphere by two rotations around X and 
H axes c-d. The initial state 1 of the system is rotated around 
the red axis first (red dotted traces) toward the intermediate 
state 2 and then rotated around the other axes (green axis 
and dashed lines) to the final state 3. The resulting combined 
rotations are shown by thick black lines. 



With the explicit demonstration of two non-commuting 
purely geometric gates, we have shown the universal char- 
acter of the realized non-adiabatic non-abelian transfor- 
mation acting on a three level quantum system. A sim- 
ilar scheme applied to two coupled three level systems 
would complete the universal set of holonomic quantum 
gates and would allow for performing all-geometric quan- 
tum algorithms, which are potentially resilient against 
noise when short pulses are useo^. Moreover, holo- 
nomic gates demonstrated for superconducting quantum 
devices, could also be applied to other three-level systems 
with similar energy level structure. 



A. Methods: Process tomography 

In order to characterize the gates, we perform full 
process tomography on the three-level system and re- 
construct the process matrix Xexp using a maximum 
likelihood procedure^. Any SU{3) operator acting on 
a three-level system can be decomposed into nine or- 
thogonal basis operators. As a full set of basis op- 
erators we choose {/qi, crg^i, -i^oi' ^01 ^ ^Oe^ ~i^Oe' 
crfg, — icr^g,£^}, where aij are Pauli operators acting 



on the levels i and j, /qi = |0)(0| + |1)(1| and 
E = |e)(e|. The process is fully determined by its 
action on a complete set of nine input states {|0), 



|e), |1), (|0) + |e))/./2, (|0) + i|e))/^/2, (|0) + |1))/V2, 
(|0)+i|l))/V2, (|e) + |l))/A (|e)+i|l))/\/2}. These 
states are prepared by sequentially applying identity, tt— , 
and 7r/2— operations (/, Rx{7t), Rx{7t/2), Ry{7r/2)) on the 
|0) ^ |e) and |e) ^ |1) transitions. After applying the 
geometric operation to the input states, we perform full 
state tomography on the respective output states —. The 
length of a typical sequence is 208 ns (five 40 ns pulses 
with 2ns separation). We calibrate the 7r-pulses on the 



|0) ^ |e) and the |e) ^ |1) transitions by measuring Rabi 
oscillations between the corresponding states. From the 
recorded 9^ = 81 measurements, the process matrix Xexp 
is reconstructed. The process is compared to the ideal one 
([2]) by calculating its fidelity as F = tr(xexpXth)- Since 
the |e) state serves only as an auxiliary state, we present 
only the reduced density matrix x, which describes the 
processes involving the |0) and |1) states and omits any 
operators acting on the |e) state. The set of basis opera- 
tors is thus given by {/, X, F, Z} = {/qi, ag^i, -icrg^, agi}. 
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